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Abstract 

The decay rates of QQ mesons {Q e c, 6) are studied in the NRQCD formahsm in terms of 
their short distance and long distance coefficients. The long distance coefficients are obtained 

\ through phenomenological potential model description of the mesons. The model parameters 

that reproduces the mass spectrum of the cc, 66 and c6 mesons are employed to study the decay 
widths of these mesons. We extract the mass spectrum as well as the reproduces the respec- 

\ tive radial wave functions from the different potential models as well as from non-relativistic 

phenomenological quark antiquark potential of the type V{r^ = — ^ + Ar^ ^ with v varying 
from 0.5 to 2. The spin hyperfine and spin-orbit interactions are employed to obtain the 

\ masses of the pseudoscalar and vector mesons. The decay constants with QCD corrections are 

computed in this model as well as in the case of other potential models for comparison. The 
digamma and dileptonic decays of cc, and 66 mesons are investigated using some of the known 
potential models without and with radiative corrections up to the lowest order. These decay 

\ width are also computed within the NRQCD formalism up to O(f^) by making uses of the 

respective spectroscopic parameters of the models. Our theoretical predictions of the decays 
of the cc, and 66 mesons and the results obtained from some of the other potential schemes 

\ are compared with the experimental values. The partial widths and life time of the Be meson 

are also computed using the model parameters and are found to be in good accordance with 
the experimental values. 

: pacs 12.39Jh, 12.40Yx, 13.20.Gd, 13.20Fc, 13.25.Hw,14.40.Nd 

. 1 Introduction 

■ Recently, there have been renewed interest in the spectroscopy of the heavy flavoured hadrons due 
to number of experimental facilities (CLEO, DELPHI, Belle, BaBar, LHCb etc) which have been 

, continuously providing and expected to provide more accurate and new informations about the 
hadrons from low flavour to heavy flavour sector [H [2] . 

The heavy flavour mesons are those in which at least one of the quark or antiquark or both 
the quark and antiquark belong to heavy flavour sector; particularly the charm or beauty. They 
are represented by QQ mesonic systems which include the quarkonia (cc and W) and Be {be or c6) 
mesons. The investigation of the properties of these mesons gives very important insight into heavy 
quark dynamics. Heavy quarkonia have a rich spectroscopy with many narrow states lying under 
the threshold of open flavour production [3],|1]. 

The success of theoretical model predictions with experiments can provide important infor- 
mation about the quark-antiquark interactions. Such information is of great interest, as it is not 
possible to obtain the QQ potential starting from the basic principle of the quantum chromody- 
namics (QCD) at the hadronic scale. In this scale it is necessary to account for non-perturbative 
effects connected with complicated structure of QCD vacuum. All this lead to a theoretical un- 
certainty in the QQ potential at large and intermediate distances. It is just in this region of large 



and intermediate distances that most of the basic hadron resonances are formed. Among many 
theoretical attempts or approaches to explain the hadron properties based on its quark structure 
very few were successful in predicting the hadronic properties starting from mass spectra to decay 
widths. For the mass predictions, the nonrelativistic potential models with Buchmiiller and Tye 
[5], Martin [6l[7l[8], Log [91110], Cornell [Tlj etc., were successful at the heavy flavour sectors while 
the Bethe-Salpeter approach under harmonic confinement [12] was successful at low flavour sec- 
tor. There exist relativistic approaches for the study of the different hadronic properties [131 [II]- 
The non-relativistic potential model has been successful for and T families, while the relativistic 
approaches yield better results in the lighter sector. Some potential models have also predicted 
the masses and various decays of the heavy-heavy mesons which are in fair agreement with the 
experimental results [151 [HI [HI [HI [13 [201 [HI [221 [231 [21] • A comprehensive review of developments 
in heavy quarkonium physics is available in ref [25] . 

The new role of the heavy flavour studies as the testing ground for the non-perturbative 
aspects of QCD, demands extension of earlier phenomenological potential model studies on quarko- 
nium masses to their predictions of decay widths with the non-perturbative approaches like NRQCD. 

The decay rates of the heavy-quarkonium states into photons and pairs of leptons are among 
the earliest applications of perturbative quantum chromodynamics (QCD) [26l|2Z]- In these anal- 
ysis, it was assumed that the decay rates of the meson factored into a short- distance part that is 
related to the annihilation rate of the heavy quark and antiquark, and long-distance factor con- 
taining all nonperturbative effects of the QCD. The short-distance factor calculated in terms of the 
running coupling constant asimq) of QCD, evaluated at the scale of the heavy-quark mass mg, 
while the long-distance factor was expressed in terms of the meson's nonrelativistic wave function, 
or its derivatives, evaluated at origin. In case of S-wave decays [2S1 [211 [2Sl [30] and in case of P- 
wave decays into photons [32], the factorization assumption was supported by explicit calculations 
at next-to- leading order in as- However, no general argument advanced for its validity in higher 
orders of perturbation theory. These divergence cast a shadow over applications of perturbative 
QCD to the calculation of annihilation rates of the heavy quarkonium states. 

In this context, an elegant effort was provided by the NRQCD formalism [33]- It consists of 
a nonrelativistic Schrodinger field theory for the heavy quark and antiquark that is coupled to the 
usual relativistic field theory for light quarks and gluons. NRQCD not only organize calculation 
of all orders in a^, but also elaborate systematically the relativistic corrections to the conventional 
formula. Furthermore, it also provides nonperturbative definitions of the long-distance factors in 
terms of matrix elements of NRQCD, making it possible to evaluate them in the numerical lattice 
calculations. Analyzing S'-wave decays within this frame work, it recover, at leading order in v'^, 
standard factorization formulae, which contain a single nonperturbative parameter. At next to 
leading order in v"^, the decay rates satisfy a more general factorization formula, which contain two 
additional independent nonperturbative matrix elements related to their radial wave functions. 

Our attempt in this paper would be then to study the heavy-heavy flavour mesons in the 
charm and beauty sector in a general frame work of the potential models. The model parameters 
used for the predictions of the masses and their radial wave functions would be used for the study 
of their decay properties using NRQCD formalism. 

For completion, we present a detail analysis of mass spectra of cc, cb and bb mesons in the 
potential scheme of coulomb plus power potential {CPP^) with the power index [u), varying from 
0.5 to 2. Spin hyperfine and spin orbit interactions are introduced to get the S-wave and P-wave 
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masses of the pseudoscalar and vector mesons. We present details of the non-relativistic treatment 
of the heavy quarks along with the computed results in section-2. The decay constants fpy of these 
mesons incorporating QCD corrections up to 0{as) are presented in section-3. The weak decay of 
the Be meson and its life time is computed in section-4, while in section-5 we present the details of 
the computations of the di-gamma decays of pseudoscalar states and the leptonic decay widths of 
the vector states of the cc and bb quarkonia in the frame work of the NRQCD formalism as well as 
other treatments incorporating different correction terms to the respective decay widths. Though 
the NRQCD formalism takes advantage of the fact that heavy quark mass is much larger than the 
other energy scales such as the binding energy scale, Aqcd and |"p^|, the energy fluctuations of 
the heavy quarks of the order of the light energy scale are implemented in pNRQCD [311 EHl 136] . 
A comprehensive comparison of the results are presented in this section. Finally we draw our 
conclusions in section-6 of this paper. 



2 Nonrelativistic Treatment for QQ systems 

Even though there are attempts based on the relativistic theory like the light front approach for 
the study of the heavy flavoured quarks, under non-relativistic approximations, they reproduces the 
results of the non-relativistic quark-potential models [37] . In the center of mass frame of the heavy 
quark- antiquark system, the momenta of quark and antiquark are dominated by their rest mass 
"^Q,Q ^ Aqcd ~ |p Ij which constitutes the basis of the non-relativistic treatment. In NRQCD, 
the velocity of heavy quark is chosen as the expansion parameter [38] . 

Hence, for the study of heavy-heavy bound state systems such as cc, cb, bb, we consider a 
nonrelativistic Hamiltonian given by [211 ESI 123] 

2 

H = M + ^ + V{r) (1) 
2m 

where 

M = tuq + iriQ, and m = — (2) 

"^Q + "^Q 

mq and niq are the mass parameters of quark and antiquark respectively, p is the relative momentum 
of each quark and V{r) is the quark antiquark potential. Though linear plus coulomb potential is 
a successful well studied non-relativistic model for heavy flavour sector, their predictions for decay 
widths are not satisfactory owing to the improper value of the radial wave function at the origin 
compared to other models [22j. Recently, we have considered a general power potential with colour 
coulomb term of the form 

V{r) = - + Ar' (3) 
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as the static quark-antiquark interaction potential {CPPy). Here, for the study of mesons, ac = '^^s, 

as being the strong running coupling constant, A is the potential parameter and z/ is a general power, 
such that the choice, u = 1 corresponds to the coulomb plus linear potential. This potential belong 
to the special choices of the generality of the potentials, V{r) = —Cr^ + Dr^ + Vq [391 SOI E] with 
Vq = a = —1, (3 = V. Choices of the power index in the range 0.5 < v < 2.0 have been explored 
for the present study. The different choices of v here, correspond to different potential forms. Thus, 
the potential parameter A can also be different numerically and dimensionally for each choices of v. 
The properties of the light-heavy flavour mesons have been calculated using the Gaussian trial wave 
function [21]. Masses and decay constant of the light-heavy systems are found to be in agreement 
with the experimental results for the choice oivK. 0.5. However, in the case of heavy-heavy systems 
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the predictions of the masses were satisfactory but the decay constants and decay rates were not 
predicted satisfactorily [22] . Hence for the present study of heavy- heavy flavour mesons, we employ 
the exponential trial wave function of the hydrogenic type to generate the Schrddinger mass spectra. 
Within the Ritz variational scheme using the trial radial wave function we obtain the expectation 
values of the Hamiltonian as {{H) = E{fi, v) ) 

E(,,.)^M^^y-[-,a.^A^±p-) (4) 

Eqn(jl]) gives the spin average mass of the ground state. For excited states the trial wave function 
is multiplied by an appropriate orthogonal polynomial function such that the excited trial wave 
function gets orthonormalized. So, it is straight forward to assume the trial wave function for the 
(n, /) state to be the form given by the hydrogenic radial wave function, 

'"'" V" f«r)' e--/' 





^■W= r2;(„ + 0! J il^ r)'e-^"^L':»Mr) (5) 

Here, /i is the variational parameter and L^^l_^{^r) is Laguerre polynomial. For a chosen value of 
V, the variational parameter, /x is determined for each state using the virial theorem 

(6) 

As the interaction potential assumed here does not contain the spin dependent part, Eqn 
(Hj) gives the spin average masses of the system in terms of the power index v. The spin average 
mass for the ground state is computed for the values of v from 0.5 to 2. We have taken the quark 
mass parameters = 4.66 GeV and = 1.31 GeV . The potential parameter A are fixed for 
each choices of v so as to get the experimental ground state masses of cc, cb and hb mesons . The 
parameters and the fitted values of A for different systems are listed in Table [TJ The experimental 
spin average masses are computed from the experimental masses of the pseudoscalar and vector 
mesons using the relation, 

MsA = Mp + -^{Mv-Mp) (7) 

For the nJ state, we compute the spin-average or the center of weight mass from the respective 
experimental values as 

_ S.2(2J+1)M., 
^^^'"-^ - E.2(2J+1) 

The fitted value of A for each case of the power index u along with other model parameters are 
tabulated in Table[T]for cc, cb and bb systems. The ground state center of weight mass of 3.068 GeV , 
6.320 GeV and 9.453 GeV are used to fit the A values for cc, cb and bb systems respectively. The 
values of A{u) thus obtained for each case of mesonic systems are then used to predict the higher S 
and P-wave masses (See Tables [2J17]) . The fitted values of A for each v in the case of the heavy- heavy 
flavour mesons are plotted in Fig ??. It is interesting to note that all the three plots intersect each 
other at v equal to 1.1 at the value of the parameter A around 0.151 GeV^"^^. It can also be seen 
that the parametric values of A for cc and cb systems are close to each other, while in the case of 
66, they are distinctly different except at z/ = 1.1. It reflects the fact that potential parameter A 
becomes independent of the distinct energy scales of these heavy mesons at around v = 1.1. 

In the case of cc and cb systems, the values of the parameter A are numerically very close to 
each other in the range of potential index 0.9 to 1.3. The predicted masses are also found to be 
in good agreement with the existing experimental states in range of power index 0.9 to 1.3 of the 
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Table 1: Parameter A for mesons in (GeV^^^) 



V 


A {cc) 


A {be) 


A {bb) 


0.5 


0.276 


0.250 


0.195 


0.7 


0.225 


0.212 


0.179 


0.9 


0.185 


0.179 


0.165 


1.0 


0.167 


0.165 


0.158 


i.i 


U.iOi 


U.ioz 


U.iOi 


1.3 


0.124 


0.129 


0.138 


1.5 


0.101 


0.109 


0.126 


1.7 


0.082 


0.092 


0.114 


1.9 


0.067 


0.077 


0.103 


2.0 


0.060 


0.070 


0.098 



ac{cc) = 0.40, ac{bc) = 0.34, ac{bb) = 0.30, 
rric = 1.31 GeV and nib = 4.66 GeV 



potential. Fig ?? shows the behaviour of |-Ri5(0)| with the potential index u for all the three (cc, 
cb and bb) mesons. Like other potential model predictions of the wave functions (at the origin) of 
be system he in between those of cc and bb systems. We obtained a model independent relationship 
similar to the one given by [12] as 

\i^,-cf ^ \A-cf^'-'^ (9) 

with q = 0.3. This relation provides the IS wave function at the origin within 2% variation for the 
choices of the potential range 0.5 < ly < 2. For 25* and 35* states we find the relation hold within 
5% for all values for u studied here. It is to be noted here that such a scaling law with smaller 
percentage variations exist here even though the potential contains coulomb part. 



2.1 Spin-Hyperfine and Spin-Orbit Splitting in Heavy-Heavy Flavour 
Mesons 

In general, the quark-antiquark bound states are represented by identified with the J 

values, with J = L + S, S = Sq + Sq, parity P = (—1)-'""'"^ and the charge conjugation C = (—1)-'""'"'^ 
and (n, L) is the radial quantum numbers. So the S'-wave {L = 0) bound states are represented by 
J^*^ = and 1 respectively. And the P-wave {L = 1) states are represented by J^'" = 1^~ 
with L = 1 and S = while J^*^ = 0++, 1++ and 2++ correspond to L = 1 and S = 1 respectively. 
Accordingly, the spin-spin interaction among the constituent quarks provides the mass splitting 
of J = O^"*" and 1 states, while the spin-orbit interaction provides the mass splitting of J^*" = 
0^^, 1^^ and 2~^~^ states. The J^'-' = state with L = 1 and S = represents the center 
of weight mass of the P-state as its spin-orbit contribution becomes zero, while the two J = 1^~ 
singlet and the J = 1^^ of the triplet P-states form a mixed sate. We add separately the spin- 
dependent part of the usual one gluon exchange potential (OGEP) between the quark antiquark for 
computing the hyperfine and spin-orbit shifting of the low-lying S'-states and P-states. Accordingly, 
the spin-spin and spin-orbit interactions are taken as ^45j 

8 » ^ ^ 

Vs^ . sJr) = Sq ■ Sq A7r5{r) (10) 

^ drnQTUQ ^ ^ ' 

... X 4:as L ■ S 

Vl ■ sir) = ^ (11) 



3 rnQfriQ r 
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Table 2: Wave function at the origin (|i?(0)|)and spin average masses of S-wave cc meson 



btate V 


/i 


1 r\\ 1 

|i?(0)| 




Exp. 


ml 

1 neory 




GeV 




(GeV) 


(GeV) 


(GeV) 


0.5 


1.068 


0.781 


3.068 






0.7 


1.177 


0.903 


3.068 






0.9 


1.264 


1.005 


3.068 






IS 1.0 


1.300 


1.049 


3.068 






1.1 


1.335 


1.090 


3.068 


3.068 


3.068" 


1.3 


1.394 


1.164 


3.068 




3.068^ 


1.5 


1.445 


1.228 


3.068 






1.7 


1.489 


1.285 


3.068 






1.9 


1.528 


1.336 


3.068 






2.0 


1.546 


1.360 


3.068 






0.5 


1.057 


0.384 


3.368 






0.7 


1.242 


0.489 


3.454 






0.9 


1.403 


0.588 


3.534 






2S 1.0 


1.473 


0.632 


3.567 






1.1 


1.540 


0.676 


3.601 


3.663 


3.662" 


1.3 


1.660 


0.757 


3.661 




3.674^ 


1.5 


1.765 


0.829 


3.713 






1.7 


1.856 


0.894 


3.756 






1.9 


1.939 


0.955 


3.796 






2.0 


1.976 


0.982 


3.814 






0.5 


1.097 


0.271 


3.550 






0.7 


1.333 


0.363 


3.712 






0.9 


1.545 


0.453 


3.870 






3S 1.0 


1.640 


0.495 


3.940 






1.1 


1.732 


0.537 


4.012 


4.040 


4.064" 


1.3 


1.901 


0.618 


4.146 




4.073^ 


1.5 


2.051 


0.692 


4.266 






1.7 


2.186 


0.762 


4.373 






1.9 


2.309 


0.827 


4.473 






2.0 


2.365 


0.857 


4.518 







Ebert(a) ^ [HI, Pandya (b) ^ 119] 
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Table 3: Derivative of wave function at the origin (|i?'(0)|) and P-wave masses of cc meson 



State 


V 


/i 


|i?'(0)| 


E[ji) 


Exp. 


Theory 






GeV 




(GeV) 


(GeV) 


(GeV) 




0.5 


1.024 


0.217 


3.313 








0.7 


1.180 


0.309 


3.373 








0.9 


1.331 


0.417 


3.426 






IP 


1.0 


1.392 


0.467 


3.450 








1.1 


1.450 


0.517 


3.473 


3.525 


3.526^^ 




1.3 


1.553 


0.614 


3.513 




3.497'' 




1.5 


1.642 


0.705 


3.547 








1.7 


1.720 


0.792 


3.576 








1.9 


1.790 


0.875 


3.603 








2.0 


1.823 


0.916 


3.615 








0.5 


1.081 


0.405 


3.519 








0.7 


1.307 


0.651 


3.666 








0.9 


1.509 


0.932 


3.808 






2P 


1.0 


1.560 


1.013 


3.872 








1.1 


1.687 


1.232 


3.936 




3.945" 




1.3 


1.847 


1.545 


4.055 




3.907'' 




1.5 


1.990 


1.862 


4.162 








1.7 


2.117 


2.174 


4.257 








1.9 


2.232 


2.481 


4.345 








2.0 


2.285 


2.631 


4.385 







Ebert(a) ^ [HI, Pandya (b) ^ [19] 
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Table 4: Wave function at the origin (|i?(0)|)and spin average masses of S-wave be meson 



btate u 


/i 


1 7~> 1 

\R{0)\ 




ml 

i neory 




GeV 




(GeV) 


(GeV) 


0.5 


1.281 


1.025 


6.320 




0.7 


1.404 


1.176 


6.320 




0.9 


1.502 


1.301 


6.320 




IS 1.0 


1.544 


1.357 


6.320 


6.317^^ 


1.1 


1.583 


1.408 


6.320 


6.319^= 


1.3 


1.652 


1.502 


6.320 




1.5 


1.711 


1.583 


6.320 




1.7 


1.763 


1.655 


6.320 




1.9 


1.807 


1.718 


6.320 




2.0 


1.825 


1.7 AA 


6.320 




0.5 


1.236 


0.486 


6.589 




0.7 


1.453 


0.619 


6.665 




0.9 


1.635 


0.739 


6.730 




2S 1.0 


1.719 


0.797 


6.761 




1.1 


1.796 


0.851 


6.790 




1.3 


1.938 


0.954 


6.844 


6.869" 


1.5 


2.061 


1.046 


6.890 


6.888^= 


1.7 


2.172 


1.132 


6.931 




1.9 


2.266 


1.206 


6.965 




2.0 


2.307 


1.239 


6.977 




0.5 


1.274 


0.339 


6.746 




0.7 


1.550 


0.455 


6.889 




0.9 


1.792 


0.565 


7.021 




3S 1.0 


1.905 


0.620 


7.085 




1.1 


2.012 


0.674 


7.147 




1.3 


2.211 


0.775 


7.265 


7.224^^ 


1.5 


2.389 


0.870 


7.372 


7.27r 


1.7 


2.550 


0.960 


7.471 




1.9 


2.692 


1.041 


7.555 




2.0 


2.755 


1.078 


7.590 





"Ebert [H]; '^Eichten [13] 
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Table 5: Derivative of wave function at the origin (|i?'(0)|) and P-wave masses of cb meson 



State z/ 




\R'm 




Theory 




GeV 




(GeV) 


(GeV) 


0.5 


1.198 


0.321 


6.542 




0.7 


1.392 


0.467 


6.597 




0.9 


1.553 


0.613 


6.641 




IP 1.0 


1.625 


0.687 


6.662 




1.1 


1.692 


0.760 


6.682 




1.3 


1.814 


0.905 


6.718 


6.749" 


1.5 


2.920 


1.042 


6.749 


6.736^= 


1.7 


2.013 


1.173 


6.777 




1.9 


2.094 


1.296 


6.799 




2.0 


2.129 


1.349 


6.806 




0.5 


1.260 


0.594 


6.720 




0.7 


1.520 


0.950 


6.851 




0.9 


1.751 


1.353 


6.969 




2P 1.0 


1.859 


1.571 


7.027 




1.1 


1.961 


1.794 


7.082 




1.3 


2.149 


2.257 


7.188 


7.145" 


1.5 


2.317 


2.725 


7.283 


7.142^= 


1.7 


2.469 


3.184 


7.370 




1.9 


2.603 


3.644 


7.445 




2.0 


2.662 


3.852 


7.476 





Ebert [H]; ^Eichten [13] 
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Table 6: Wave function at the origin (|i?(0)|)and spin average masses of S-wave bb meson 



State z/ 


/i 


\R{0)\ 




Theory 




GeV 




(GeV) 


(GeV) 


0.5 


1.985 


1.977 


9.453 




0.7 


2.122 


2.186 


9.453 




0.9 


2.238 


2.368 


9.453 




IS 1.0 


2.288 


2.447 


9.453 




1.1 


2.336 


2.525 


9.453 


9.445'^ 


1.3 


2.402 


2.662 


9.453 


9.453*^ 


1.5 


2.491 


2.780 


9.453 




1.7 


2.554 


2.885 


9.453 




1.9 


2.611 


2.982 


9.453 




2.0 


2.638 


3.030 


9.453 




0.5 


1.701 


0.784 


9.701 




0.7 


1.979 


0.984 


9.758 




0.9 


2.227 


1.175 


9.812 




2S 1.0 


2.338 


1.264 


9.838 




1.1 


2.442 


1.349 


9.861 


10.016'^ 


1.3 


2.636 


1.513 


9.905 


10.008"^ 


1.5 


2.807 


1.663 


9.944 




1.7 


2.958 


1.790 


9.977 




1.9 


3.094 


1.924 


10.008 




2.0 


3.158 


1.984 


10.023 




0.5 


1.692 


0.519 


9.826 




0.7 


2.053 


0.694 


9.935 




0.9 


2.385 


0.868 


10.043 




3S 1.0 


2.538 


0.953 


10.095 




1.1 


2.684 


1.036 


10.144 


10.348*^ 


1.3 


2.957 


1.199 


10.241 


10.351"^ 


1.5 


3.205 


1.353 


10.330 




1.7 


3.427 


1.496 


10.410 




1.9 


3.631 


1.631 


10.485 




2.0 


3.727 


1.696 


10.522 





Ebert (a) ^ [H], Gupta (d) ^ [H] 
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Table 7: Derivative of wave function at the origin (|i?'(0)|) and P-wave masses of bb meson 



State z/ 




\R'm 




Exp 


Theory 




GeV 




(GeV) 


(GeV) 


(GeV) 


0.5 


1.654 


0.718 


9.670 






0.7 


1.904 


1.021 


9.712 






0.9 


2.121 


1.337 


9.751 






IP 1.0 


2.218 


1.446 


9.768 






1.1 


2.308 


1.653 


9.784 


9.900 


9.901" 


1.3 


2.474 


1.966 


9.815 




9.900"^ 


1.5 


2.621 


2.271 


9.842 






1.7 


2.750 


2.559 


9.864 






1.9 


2.866 


2.838 


9.885 






2.0 


2.921 


2.977 


9.896 






0.5 


1.669 


1.200 


9.808 






0.7 


2.016 


1.922 


9.908 






0.9 


2.333 


2.770 


10.006 






2P 1.0 


2.478 


3.223 


10.053 






1.1 


2.617 


3.692 


10.097 


10.260 


10.261" 


1.3 


2.876 


4.677 


10.184 




10.258^^ 


1.5 


3.111 


5.691 


10.264 






1.7 


3.321 


6.699 


10.335 






1.9 


3.513 


7.708 


10.401 






2.0 


3.604 


8.218 


10.434 







Ebert (a) ^ [H], Gupta (d) ^ [H] 
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Table 8: S-Wave and P-Wave Masses (in GeV) of cc meson 



u 




1 


's. 


l^Pi 


1 ''Po 


1 ^Pi 


0.5 


3 


000 


3 


092 


3.313 


3.292 


3.302 


0.7 


2 


980 


3 


100 


3.373 


3.341 


3.357 


0.9 


2 


960 


3 


109 


3.429 


3.383 


3.406 


1.0 


2 


950 


3 


112 


3.450 


3.398 


3.424 


1.1 


2 


942 


3 


116 


3.473 


3.414 


3.444 


1.3 


2 


926 


3 


123 


3.513 


3.441 


3.477 


1.5 


2 


912 


3 


129 


3.547 


3.461 


3.504 


1.7 


2 


899 


3 


134 


3.576 


3.477 


3.526 


1.9 


2 


887 


3 


141 


3.603 


3.491 


3.547 


2.0 


2 


882 


3 


144 


3.615 


3.497 


3.556 



1 'P, 



2'Sn 



2 'Si 2'Po 2 •^Pi 2 ^P2 3^50 



3.323 
3.389 
3.451 
3.477 
3.503 
3.550 
3.590 
3.625 
3.658 
3.673 



3.352 
3.427 
3.495 
3.522 
3.549 
3.597 
3.636 
3.668 
3.696 
3.708 



3.375 
3.464 
3.547 
3.583 
3.619 
3.683 
3.739 
3.788 
3.832 
3.852 



3.519 


3.494 


3.507 


3.666 


3.623 


3.644 


3.808 


3.742 


3.775 


3.872 


3.792 


3.832 


3.936 


3.843 


3.889 


4.055 


3.933 


3.994 


4.162 


4.009 


4.085 


4.257 


4.073 


4.165 


4.345 


4.129 


4.237 


4.385 


4.153 


4.269 



3.531 
3.687 
3.842 
3.911 
3.983 
4.116 
4.239 
4.394 
4.453 
4.501 



3.541 
3.697 
3.846 
3.912 
3.979 
4.102 
4.212 
4.309 
4.396 
4.436 



mm 


2.980 


3.097 


3.511 


3.415 




3.556 


3.622[50J 


3.6^ 


^6 










3.929 




4 


m 


2.979 


3.096 


3.526 


3.424 


3.511 


3.556 


3.588 


3.6^ 


^6 


3.945 


3. 


854 


3.929 


3.972 


3.991 


4 


m 


2.980 


3.097 


3.527 


3.416 


3.508 


3.558 


3.597 


3.6^ 


^6 


3.960 


3. 


844 


3.894 


3.994 


4.014 


4 



The value of the radial wave function R(0) for 0~ and 1 states would be different due to their 
spin dependent hyperfine interaction. The spin hyperfine interaction of the heavy flavour mesons 
are small and this can cause a small shift in the value of the wave function at the origin. Thus, 
many other models do not consider this contribution to their value of R(0). However, we account 
this correction to the value of R(0) by considering 



P„j(0) = P(0) 



1 + iSF)j^ 



MsA 



(12) 



Where {SF)j and < Esd >nj is the spin factor and spin interaction energy of the meson in the 
nJ state, while R{0) and Msa correspond to the radial wave function at the zero separation and 
spin average mass respectively of the QQ system. It can be seen that Eqn (fT^ provides the average 
radial wave function given by [33] as 

P(0) = 5p±^ (13) 

It is found that the computed mass increases with increase of u. The computed results for the 
pseudoscalar(P) and Vector (V) mesons in the case of cc, cb, bb systems are tabulated in Tables 
[SlfTUl The spin-spin hyperfine and spin-orbit interactions are computed perturbatively to get the 
masses of rjc, J/ip^ Pc, Pci Vby and T states. The results are compared with known experimental 
values as well as with other theoretical predictions. Mass predictions with u between 1.0 and 1.5 
are found in accordance with the experimental results [2]. 



3 Decay constants (fp/v) of the heavy flavoured mesons 

The decay constants of mesons are important parameters in the study of leptonic or non-leptonic 
weak decay processes. The decay constants of pseudoscalar (/p) and vector (/y) mesons are obtained 
by prarameterizing the matrix elements of weak current between the corresponding mesons and the 
vacuum as 

{0\Qi''i,Q\P,ik))=tfpk'' (14) 
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Table 9: S-Wave and P-Wave Masses (in GeV) of cb meson 





i Do 


i Dl 


J- -< 1 


-L -T 


J- -T 1 


1 3r, 
-L -(2 


Z Do 


93c 

Z Dl 


91 P, 


93 Pr. 


9 3 R 
Z 11 


9 3 p„ 
Z 12 


Dq 


3 C 
Dl 




Fi 991 


6 330 


6 542 


6 534 


6 538 


6 546 


6 582 


6 591 


6 790 


6 71 1 


6 71 5 


6 796 


6 743 


6 747 


7 


fi 283 


6 334 


6 597 


6 584 


6 590 


6 603 


6 655 


6 669 


6 851 


6 834 


6 840 


6 859 


6 884 


6 891 


u.y 


OVQ 

o.z / o 


D.OOO 


o.o4i 


D.D/4 


D.boo 


O.ObU 


0. / ib 


0. ^ ob 


D.yoy 


P. C\A A 

0.944 


D.yb / 


o.yoz 


V ni 
i .\jlZ 


/ .Uz4 


1.0 


6.269 


6.337 


6.662 


6.642 


6.652 


6.672 


6.743 


6.767 


7.027 


6.997 


7.012 


7.042 


7.075 


7.089 


1.1 


6.265 


6.338 


6.682 


6.659 


6.671 


6.693 


6.770 


6.797 


7.082 


7.047 


7.065 


7.099 


7.135 


7.151 


1.3 


6.259 


6.341 


6.718 


6.691 


6.704 


6.732 


6.819 


6.852 


7.188 


7.142 


7.165 


7.211 


7.249 


7.271 


1.5 


6.252 


6.344 


6.749 


7.716 


7.733 


7.765 


6.860 


6.900 


7.283 


7.225 


7.254 


7.312 


7.351 


7.379 


1.7 


6.247 


6.347 


6.777 


7.739 


7.758 


7.796 


6.896 


6.943 


7.370 


7.300 


7.335 


7.405 


7.445 


7.479 


1.9 


6.241 


6.348 


6.799 


7.756 


7.777 


7.820 


6.924 


6.978 


7.445 


7.363 


7.404 


7.486 


7.525 


7.565 


2.0 


6.237 


6.348 


6.806 


7.762 


7.784 


7.829 


6.935 


6.991 


7.476 


7.388 


7.432 


7.519 


7.558 


7.601 


m 


6.270 


6.332 


6.749 


6.699 


6.734 


6.762 


6.835 


6.881 


7.145 


7.091 


7.126 


7.145 


7.193 


7.235 


m 


6.256 


6.337 


6.755 


6.700 


6.730 


6.747 


6.899 


6.929 


7.169 


7.108 


7.135 


7.142 


7.280 


7.308 



Table 10: S-Wave and P-Wave Masses (in GeV) of hh meson 



l^^o 1 ^^1 l^A 1 ^Po 1 ^Pi 1 ^P2 2^So 2 3^1 2iPi 23Po 2 



2 3P2 



3"^ 5*0 



0.5 
0.7 
0.9 
1.0 
1.1 
1.3 
1.5 
1.7 
1.9 
2.0 



9.426 
9.419 
9.414 
9.411 
9.408 
9.403 
9.399 
9.394 
9.390 
9.389 



9.463 
9.465 
9.467 
9.468 
9.468 
9.470 
9.472 
9.473 
9.474 
9.475 



9.672 
9.716 
9.757 
9.775 
9.791 
9.824 
9.852 
9.877 
9.900 
9.911 



9.664 
9.703 
9.740 
9.755 
9.769 
9.797 
9.820 
9.840 
9.857 
9.866 



9.670 
9.712 
9.751 
9.768 
9.784 
9.815 
9.842 
9.864 
9.885 
9.896 



9.683 
9.731 
9.774 
9.792 
9.809 
9.840 
9.866 
9.887 
9.905 
9.913 



9.696 
9.751 
9.803 
9.826 
9.846 
9.888 
9.924 
9.955 
9.982 
9.995 



9.702 

9.760 

9.816 

9.841 

9.865 

9.910 

9.951 

9.985 

10.017 

10.032 



9.808 

9.908 

10.006 

10.053 

10.097 

10.184 

10.264 

10.335 

10.401 

10.434 



9.803 

9.898 

9.991 

10.035 

10.076 

10.155 

10.228 

10.291 

10.350 

10.379 



9.806 

9.903 

9.999 

10.044 

10.086 

10.170 

10.246 

10.313 

10.376 

10.406 



9.811 

9.913 

10.014 

10.062 

10.108 

10.198 

10.282 

10.357 

10.428 

10.462 



mm 




9.460 




9.; 


860 


9.; 


893 


9.913 




10.023 




10.232 


10.255 


10.268 




m 


9.400 


9.460 


9.901 


9.; 


863 


9.; 


892 


9.913 


9.993 


10.023 


10.261 


10.234 


10.255 


10.268 


10.328 


m 


9.414 


9.461 


9.900 


9.; 


861 


9.; 


891 


9.912 


9.999 


10.023 


10.262 


10.231 


10.255 


10.272 


10.345 
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{0\QrQ\Vik,e)) = fvMve'' (15) 

where k is the meson momentum, and My are the polarization vector and mass of the 
vector meson. In the relativistic quark model, the decay constant can be expressed through the 
meson wave function ^py{p) in the momentum space as [H] 



fp,\ 



12 



M, 



py 



2Eq{p) J\\ 2Eq{p) 



' Eq{p) +m. 




1 + A 



p 



py 



[EQ{p)+mQ][EQip) 




py\ 



with Ap = —1 and Ay = —1/3. In the nonrelativistic limit 



p 



0, this expression reduces 



to the well known relation between fpy and the ground state wave function at the origin i/jpy{0), 
the Van-Royen-Weisskopf formula |l6j. Though most of the models predict the mesonic mass 
spectrum successfully, there are disagreements in the predictions of the pseudoscalar and vector 
decay constants. For example, most of the cases, the ratio was predicted to be > 1 as mp < my 
and their wave function at the origin ipp{0) ~ ipv{0) [SI- The ratio computed in the relativistic 
models [17] predicted the ratio < 1, particularly in the heavy flavours sector. The disparity 
of the predictions of these decay constants play decisive role in the precision measurements of the 
weak decay parameters. So, we re-examine the predictions of the decay constants under different 
potential schemes discussed in the present work. Incorporating a first order QCD correction factor 
, we compute them using the relation. 



fp/v 



12 



^p/y(0) 



M, 



p/v 



(17) 



Where C{as) is the QCD correction factor given by [i8l H9] 



vr 



In 

niQ + niQ m, 



Here 5^ = 2 and 5^ = 8/3. The computed fp and /y for cc, cb and bb systems using Eqn([T] 
& (ITS]) and the predicted radial wave functions at the origin i?„j(0) of the respective mesons are 
tabulated in Tables [TT] - [T^ The decay constants without and with the QCD corrections are also 
listed as fpy and fpy{cor.) in the table. The plot of fp vs {mq + mg) shows (see fig 3) deviations 
from linearity as against the predictions of a linear scaling between the weak decay constant and 
the sum of quark antiquark masses justified within a renormalized light front QCD inspired model 
for quark antiquark bound states j51j . 



4 Weak decay of meson 

The Decay properties of B'^ (he) meson is of interest as it decays only through weak interactions 
[m [211 [52], [53]. This is due to the fact that its ground state energy lies below the (BD) threshold 
and has non vanishing flavour. This eliminates the uncertainties encountered due to strong decays 
and provides a clear decay width and lifetime for 5+ meson, which helps to fix more precise value 
of the weak decay parameters such as the CKM mixing matrix elements (V^b, V^^) and the leptonic 
decay constant {fp). Adopting the spectator model for the charm-beauty system [21], the total 
decay width of B'^ meson can be approximated as the sum of the widths of 6-quark decay keeping 
c-quark as spectator, the c-quark decay with 6-quark as spectator, and the annihilation channel 
B'^ /+z/;(cs, Ms), I = e,fi,T with no interference assumed between them. 
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Table 11: Decay constants(/p & fv) (in MeV) of IS* and cc mesons states. 



Models 


RpiO) 


RviO) 


fp 


fp{cor.) 


fv 


fv{cor.) 








MeV 


MeV 


MeV 


MeV 


ERHM 


0.726 


0.752 


410 


317 


418 


323 


BT 


0.874 


0.909 


499 


382 


505 


389 


PL 


0.971 


1.009 


550 


399 


560 


407 


LOG 


0.877 


0.914 


496 


379 


506 


387 




1.171 


1.217 


663 


532 


676 


543 




763 


787 


430 


348 


437 


326 


7 


875 


912 


495 


401 


506 


377 


0.9 


0.967 


1.018 


549 


444 


564 


421 


1.0 


1.005 


1.063 


572 


463 


589 


439 


1.1 


1.041 


1.107 


593 


480 


613 


457 


1.3 


1.104 


1.184 


627 


507 


655 


488 


1.5 


1.158 


1.252 


663 


536 


692 


516 


1.7 


1.204 


1.311 


691 


559 


724 


539 


1.9 


1.245 


1.366 


716 


579 


753 


561 


2.0 


1.264 


1.391 


728 


589 


767 


571 










335± 


459± 


416± 










75 [55] 


28[56] 


6[55] 



Table 12: Pseudoscalar and Vector meson decay constants(/p & fv) (in MeV) of IS' be meson state. 





RpiO) 


Rv{0) 


fp 


fp{cor.) 


fv 


fv{cor.) 


u=Q.5 


1.021 


1.027 


398 


356 


399 


336 


0.7 


1.169 


1.178 


456 


408 


457 


385 


0.9 


1.291 


1.304 


504 


451 


506 


426 


1.0 


1.346 


1.361 


525 


470 


528 


445 


1.1 


1.396 


1.413 


545 


488 


548 


461 


1.3 


1.487 


1.507 


581 


520 


585 


492 


1.5 


1.565 


1.589 


612 


548 


616 


519 


1.7 


1.635 


1.662 


639 


572 


645 


542 


1.9 


1.695 


1.725 


663 


594 


669 


563 


2.0 


1.722 


1.758 


674 


603 


682 


574 










433 [S] 


503115 


418 
± 24 [561 
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Table 13: Decay constants(/p & fy) (in MeV) of IS bb meson state. 



Models 


RJO) 


RJO) 


fp 


fj:>(cOT 1 


fv 

J V 


f\/(cOT 1 








MeV 


MeV 


MeV 


MeV 


ERHM 


2.232 


2.235 


709 


601 


710 


601 


BT 


2.527 


2.551 


807 


683 


810 


686 


PL 


2.132 


2.146 


680 


563 


682 


565 


LOG 


2.206 


2.221 


703 


594 


706 


596 


Cornell 


3.706 


3.762 


1185 


1022 


1194 


1029 


u=0.5 


1.971 


1.979 


627 


537 


629 


509 


0.7 


2.178 


2.189 


693 


594 


695 


563 


0.9 


2.358 


2.371 


751 


643 


753 


609 


1.0 


2.436 


2.451 


776 


665 


778 


630 


1.1 


2.513 


2.529 


801 


686 


803 


650 


1.3 


2.648 


2.667 


844 


723 


847 


685 


1.5 


2.764 


2.785 


881 


755 


884 


715 


1.7 


2.867 


2.891 


914 


783 


918 


743 


1.9 


2.962 


2.989 


945 


809 


949 


768 


2.0 


3.009 


3.037 


960 


822 


964 


780 



711 [19J 708 
± 8132 



Table 14: Decay widths (in 10 ^ eV) and lifetime r (in ps) of B'^ meson 



Model 


r[Anni) 


r(5. 




r (in 


PS) 




a 


b 


a 


h 


a 


b 


u = 0.5 


0.370 


0.370 


12.47 


13.18 


0.530 


0.499 


0.7 


0.486 


0.484 


12.59 


13.30 


0.523 


0.495 


0.9 


0.596 


0.593 


12.70 


13.41 


0.518 


0.491 


1.0 


0.644 


0.642 


12.75 


13.46 


0.516 


0.489 


1.1 


0.693 


0.690 


12.79 


13.51 


0.515 


0.487 


1.3 


0.786 


0.783 


12.89 


13.60 


0.511 


0.484 


1.5 


0.871 


0.867 


12.89 


13.68 


0.507 


0.481 


1.7 


0.951 


0.950 


12.97 


13.76 


0.504 


0.478 


1.9 


1.023 


1.020 


13.05 


13.83 


0.502 


0.476 


2.0 


1.053 


1.050 


13.15 


13.87 


0.500 


0.475 


m 










o.46l^:l« 




m 




1.40 




14.00 


0.47 




[52] 




0.67 




8.8 


0.75 





16 



Accordingly, the total width is written as [21] 



r(5c -^X) = T{b ^X) + r(c X) + T{Anni) (19) 

Neglecting the quark binding effects, we obtain for the b and c inclusive widths in the spectator 
approximation as [24j 



r{b^x) = 



9 Gl\V,bfml 

1927r3 
= 7.97 X lO-^el^(a) 

= 8.66 X 10-^e1/(6) (20) 



r(c^X) 



5 Gl\V,sfml 



1927r3 
4.13 X 10-^eV{a) 

4.15 X IQ-'^evlb) (21) 



Here we have used the model quark masses and the two values (a) and (b) correspond to the two 
set of values for the CKM matrix elements (a)— >|Vcs| = 0.97296, = 0.04221 as used in reference 
[1] and (b)— \Vcs\ = 0.975, \Vcb\ = 0.044 as the upper bound provided by particle data group. The 
values of r{B X) and r(c — » X) in Bethe-Salpeter model [24j and relativized quark model [52] 
are 7.5 & 5.1 and 4.8 & 3.3 (widths are in 10""^ eV) respectively. 

Employing the computed mass of the I^Sq state {Mb J and /^^ values obtained from the 
present study, the width of the annihilation channel is computed using the expression given by [24], 

V{Anm) = ^ \V,f fl^M,,m] [l - ^ j C„ (22) 

2 

Where = 3 \ Vcs\ for cs, and niq is the mass of the heaviest fermions. 

The computed widths and lifetime in our CPPy model are listed in Table [TH Our predictions 
for the life time with the potential index 0.5 ^ z/ ^ 2 lie well within the experimental error bar. 



5 Decay rates of quarkonia 

Along with the mass spectrum, successful predictions of various decay widths of heavy flavoured 
systems have remained as testing ground for the success of phenomenological models. Experimen- 
tally, the excited states and the leptonic, di-gamma and other hadronic decay width, of the heavy 
flavour mesons have been reported. However, experimentally, the pseudoscalar bb bound state rji, 
is still elusive though experimental search for this state at the di-gamma decay channel has been 
initiated recently 07]. 

As an attempt to improve the theoretical predictions involving the phenomenological descrip- 
tion of the meson, using the redial wave functions and other model parameters of the different 
potential models we study the decay of ^5*0 quarkonium into di-gamma and the decay of ^5*1 into 
lepton pairs using the NRQCD formalism [SB]- It is expected that the NRQCD formalism has 
all the corrective contributions for the right predictions of the decay rates. NRQCD factorization 
expressions for the decay rates of quarkonium and decay are given by [31] 
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rCSo ^ 77) 



rrir 



<oixvr^o> 



2 Gjj(^^Sq 



-Re 



m 



Q 



my 



So\^ph\o><o\x\-^'^ri^\'So> 



(23) 



T^Si ^ e+e-) 



Fee^^i) 



rl /I ^ 



2 Gee('^l) „ 

+ \ Re 



S^\^P^ax\0 X 0\x^a{--'Sfi,\'S^ > 



Heei So) 



Re 



<3 S^\^^ax\0 X 0\x^ai--'^)^^pfS, > 



(24) 



The short distance coefficients F's and G's of the order of and are given by [5Tj 

^2 \ 



-a 



1367r 
45 



(25) 

(26) 
(27) 



-Fee( 'S'l) 



+ 



,a 



3 - TT 

-117.46 + 0.82n^ + ^^/n( — ; 

27 /iA 



71 



(28) 



G^ee(^'S'i) 



87rQ2 2 



-a 



2„ 2 



(29) 
(30) 



The matrix elements that contributes to the decay rates of the S wave states into riQ — >■ 77 and 
ip — > e~^e~ through next-to- leading order in v"^, the vacuum-saturation approximation gives [33] 



(31) 



(32) 



<^ 5o|Pi('5o)r5o >= Re[<' So\iJ^x\0 > 
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<0\xH-'^'^ri^\'So>]+0{v'T) (33) 



S,\V,i'S,)\^S, >= Re[<' S^\iJ^ax\0 > 

<0\x^ai-'-'nM'S,>]+O{v'T) (34) 



So\Q\CSo)\'So >=< 0\xK-'-'^yi^\'So > (35) 



Si\QieSi)fS^ >=< 0\x^a{-'-'SM'S^ > (36) 



So\QlCSo)\'So >=< 0\x\-'-^)'i^\'So > (37) 



Si\QieS,)fS^ >=< 0|xV(-^:^)Vl'^i > (38) 

The Vacuum saturation allows the matrix elements of some four fermion operators to be expressed 
in terms of the regularized wave-function parameters given by 



So\OCSo)\'So>=^\Rpm' (39) 
S,\OCS,)\'S^>=^\Rvm' (40) 



So\ViCSo)\'So >= -—\R*p v'Rp\ (41) 

Ztt 



<3 S,\V,eS,)\'S, >= -—\R*v V'Rvl (42) 

Zn 



<' So\QlCSo)\'So >= -v/|^V2i?p (43) 



<^ S^\QlCS^)\'S^ >= -jAv^i^v- (44) 



So\QlCSo)\'So >= ^V'iV^p) (45) 



S,\QlCS,)\'S, >= ^V'iV'Ry) (46) 

ZTT 
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Table 15: Decay rates (in keV) of +-^77 and the relevant correction terms of rjc and rib mesons. 



Systems Models 


To 


Tr 


F 


^NRQCD 


^NRQCDfrs 


others 










iinto D( 


nnto 






ERHM 


7.460 


-2.855 


4.605 


4 005 


4 225 


— 




BT 


10.870 


-4.206 


6.664 


6.555 


6.561 


— 


7.2±0.7±2.0[1] 


PL (Martin) 


13.406 


-6.196 


7.210 


8.434 


10.691 


— 


7.500P1 


Vc Log 


10.937 


-4.349 


6.588 


6.691 


6.697 


— 




Cornell 


19.512 


-6.581 


12.931 


13.779 


17.447 


— 


9.02 ± 0.8 [33] 


CPP^ =0.5 


8.173 


-2.635 


5.538 


2.511 


6.078 


2.992 




0.7 


10.918 


-3.521 


7.397 


3.706 


7.810 


4.087 




0.9 


13.465 


-4.342 


9.123 


4.925 


9.391 


5.102 




1.0 


14.649 


-4.724 


9.925 


5.552 


10.077 


5.549 




1.1 


15.812 


-5.099 


10.713 


6.171 


10.761 


6.012 




1.3 


17.987 


-5.800 


12.187 


7.387 


12.016 


7.095 




1.5 


19.971 


-6.440 


13.531 


8.556 


13.142 


7.536 




1.7 


22.788 


-7.026 


15.762 


9.700 


14.166 


8.170 




1.9 


23.502 


-7.578 


16.924 


10.789 


15.139 


8.736 




2.0 


24.297 


-7.835 


16.462 


11.295 


15.596 


9.006 




ERHM 


0.444 


-0.114 


0.326 


0.315 


0.317 


— 




BT 


0.574 


-0.149 


0.424 


0.445 


0.455 


— 




PL (Martin) 


0.406 


-0.118 


0.288 


0.312 


0.340 


- 


0.364 [33] 


Vb Log 


0.435 


-0.115 


0.320 


0.337 


0.345 




0.490 [59] 


Cornell 


1.244 


-0.290 


0.954 


1.015 


1.112 






CPP^ =0.5 


0.345 


-0.086 


0.259 


0.254 


0.254 


0.195 




0.7 


0.422 


-0.106 


0.316 


0.310 


0.311 


0.256 




0.9 


0.495 


-0.124 


0.371 


0.365 


0.366 


0.321 




1.0 


0.529 


-0.132 


0.397 


0.390 


0.391 


0.353 




1.1 


0.563 


-0.141 


0.422 


0.416 


0.416 


0.386 




1.3 


0.626 


-0.156 


0.470 


0.462 


0.463 


0.435 




1.5 


0.683 


-0.171 


0.512 


0.505 


0.505 


0.510 




1.7 


0.735 


-0.184 


0.551 


0.544 


0.545 


0.570 




1.9 


0.786 


-0.196 


0.590 


0.582 


0.582 


0.628 




2.0 


0.811 


-0.203 


0.608 


0.600 


0.601 


0.657 





ERHM [T9l[20j, BT ^, PL (Martin) [6J, Log [9], Cornell ^ 
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Table 16: Decay rates (in keV) of 1 ^ I and the relevant correction terms of J/ip and T 
mesons. 



bystems Models 


i vw 


rad 


i 


'^NRQCD 


i NRQCDfrs 












upto O(f^) 


upto 0[v^) 






ERHM 


5.595 


-3.381 


2.214 


2.543 


3.246 






J/ip BT 


8.152 


-4.982 


3.170 


2.539 


2.809 




5.55±0.14± 0.02 


rL (iVlartmj 


10.055 


O /I "1 

-7.341 


2.714 


3.311 


4.698 






Log 


8.203 


-0.171 


3.057 


1.967 


2.094 






Cornell 


14.634 


-7.701 


6.933 


7.920 


10.294 






CPPy= 0.5 


6.130 


-0.624 


5.506 


4.212 


4.973 


0.973 




0.7 


8.189 


-0.845 


7.344 


6.199 


7.701 


1.676 




0.9 


10.153 


-1.065 


9.088 


8.320 


10.815 


2.447 




1.0 


11.053 


-1.165 


9.888 


9.353 


12.398 


2.822 




1.1 


11.946 


-1.268 


10.678 


10.430 


14.089 


3.227 




1.3 


13.621 


-1.463 


12.158 


12.558 


17.550 


3.996 




1.5 


15.165 


-1.645 


13.520 


12.605 


19.037 


4.749 




1.7 


16.582 


-1.813 


14.769 


16.643 


24.587 


5.467 




1.9 


17.920 


-1.982 


15.938 


18.659 


28.232 


6.185 




2.0 


18.549 


-2.061 


16.488 


19.634 


30.032 


6.534 




ERHM 


1.320 


-0.540 


1.303 


1.221 


1.228 






T B.T. 


1.720 


-0.076 


1.644 


1.249 


1.267 






rL(Martmj 


1.218 


-0.761 


0.457 


0.693 


0.774 




1.340 


Log 


1.305 


-0.032 


1.273 


0.924 


0.943 




±0.018 


Cornell 


3.733 


-0.232 


3.501 


2.025 


2.270 






CPP^= 0.5 


1.035 


-0.010 


1.025 


0.935 


0.938 


0.710 


1.43[SI] 


0.7 


1.266 


-0.013 


1.253 


1.144 


1.148 


0.933 




0.9 


1.485 


-0.015 


1.470 


1.344 


1.349 


1.165 




1.0 


1.587 


-0.017 


1.570 


1.436 


1.442 


1.279 




1.1 


1.690 


-0.018 


1.678 


1.529 


1.535 


1.397 




1.3 


1.878 


-0.020 


1.858 


1.702 


1.709 


1.575 




1.5 


2.047 


-0.022 


2.025 


1.857 


1.865 


1.844 




1.7 


2.206 


-0.024 


2.182 


2.002 


2.010 


2.057 




1.9 


2.357 


-0.026 


2.331 


2.141 


2.149 


2.266 




2.0 


2.433 


-0.026 


2.407 


2.210 


2.219 


2.371 





ERHM [iniEO], BT [5], PL (Martin) [6], Log [9J, Cornell [11] 
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9 

We have computed the V Rp/^ term as per ref [57]. Accordingly, 

V^i? = es i?y , as r (47) 

where is the binding energy and M is the mass of the respective mesonic state. The binding 
energy is computed as = M — (2mQ). The RHS of the Eqn (H5l) and (H6l) are computed by 
assuming that < >^~< >. For comparison, we also compute the decay widths with the 
conventional V-W formula with and without the radiative corrections. 

Accordingly, the two photon decay width of the pseudoscalar meson is given by [22j 

r(o-+^27) = ro + rR (48) 

Here Fq is the conventional Van Royen-Weisskopf term for the ^ 77 decays [16], where 
is due to the radiative corrections for this decay which is given by 



To = RUO) (49) 



and 

(50) 




Similarly, the leptonic decay width of the vector meson is computed as 

Y(r--^l+l-)=TvW + Trad (51) 

where 



vw 



1^ d2 



KiO) (52) 



Trad, the radiative correction is given by 



rad 



-^as Tvw (53) 

OTT 

It is obvious to note that the computations of the decay rates and the radiative correction 
term described here require the right description of the meson state through its radial wave function 
at the origin, R(0) and its mass M along with other model parameters like and the model quark 
masses. Generally, due to lack of exact solutions for colour dynamics, Rp/y(0) and M are considered 
as free parameters of the theory [5Zj. However, it is appropriate to employ the phenomenological 
model spectroscopic parameters such as of the predicted mesonic mass and the corresponding wave 
function for the computations of the decay widths. In many cases of potential model predictions, 
the radial wave function at the origin are over estimated as for the decay rates are concerned. In 
such cases, it is argued that the decay of QQ occurs not at zero separation, but at some finite 
Q — Q radial separation. Then arbitrary scaling of the radial wave function at zero separation are 
done to estimate the decay rates correctly [ll] . In the present computation of the decay rates using 
the NRQCD formalism we present our results obtained by using the radial wave function and their 
derivatives at zero separation {T]^rqcd) as well as at a finite radial separation of r,,, (Fjvrqcd/^J- 
We defined by 

Mpiv 
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of the mesonic state. It is similar to the compton radius and we call it as color compton radius of 
the QQ systems. Here, A^c = 3 and cq is the charge of the quark in terms of the electron charge. 

The computed decay widths for 0^^ ^77, are presented in Table fTSl and for 1 are 
listed in Table fT6l In the case of Tj^rqcd terms up to 0(i;^) and terms up to 0{v'^) are separately 
tabulated to highlight their contributions in the respective decays. 

6 Conclusion and discussion 

In this paper, we have made a comprehensive study of the heavy-heavy flavour mesonic systems 
in the general frame work of potential models. The potential model parameters and the masses of 
the charmed and beauty quark obtained from the respective quarkonia mass predictions have been 
employed to study their decay properties in the frame work of NRQCD formalism as well as using 
the conventional Van-Royen-Weisskopf nonrelativistic formula. We have also made parameter free 
prediction of the weak decay properties of meson. The weak decay constants of the pseudoscalar 
(/p) and the vector meson {fy) computed here are is found to be in accordance with the recent 
predictions based on relativistic Bethe- Salpeter method [56]. The departure from the predicted 
linear dependence of fp with the mesonic masses within the effective light-front model in the heavy 
flavour sector suggest the requirement of more refined mechanism related to their wavefunctions 
incorporating the confinement and hyperfine splitting. 

Masses of the pseudoscalar and vector mesons and the values of the radial wave function at 
the origin for cc, cb and hh systems are computed in different potential schemes. The respective 
decay constants (/p, fv) are computed with and without QCD corrections. Using the predicted 
masses and redial wave functions at the origin, the digamma, leptonic, light hadronic decays of 
quarkonia and the weak decay properties of B'^ mesons are studied. For the mass predictions and 
for the decay rates the present results based on {CPPu) are found to be in accordance with other 
potential model predictions as well as with the experimental values. 

The theoretical [CPPy] predictions of the decay widths for J/ip ^ l^l^ and T — >■ /+/^ as 
presented in Table [16] are found to be in accordance with other potential model predictions with 
the radiative correction as well as with the widths computed using NRQCD formalism. 
Though the radiative corrections are found to be important in most of the phenomenological mod- 
els, the NRQCD predictions with their matrix elements computed at finite radial separation defined 
through the 'color compton radius' are found to be in better agreement with the experimental values 
for most of the cases. 

It is interesting to note that the ERHM[20] predictions of the di-gamma decay widths of rjc 
and leptonic decay widths of J/ip and T are in good agreement with the respective experimental 
results with out any correction to the Van-Royen-Weisskopf formula. 

The NRQCD width for rjc —>■ 77 predicted in the present study based on the potential model 
parameters of BT [5], Log [9], CPPi.=0.7,0.9 are close to the experimental value of 7.2±0.7±2.0 
keV reported by PDG2006 [1]. However for the rjf, — > 77 case, most of the model predictions based 
on NRQCD formalism are very close to similar theocratical predictions of [33]. The predictions 
based on V-W formula with radiative corrections are also found to be in close agreement with the 
prediction of [33] and [59] respectively. 
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The predictions of rjb mass spectra, its hyperfine mass split (T — rib) of 60 Mel^, its decay 
constant fp and the digamma width etc are important for the experimental hunting of 77^ state. 

In the case of the dileptonic width of cc state, our predictions based on the NRQCD formalism 
with the finite range correction for the inter quark potential index 1.5 < z/ < 1.7 are in fair agrement 
with the experimental value of 5.55±0.14 keV; while that bb system the NRQCD prediction is 
in good agrement with value of 1.340 ± 0.018 keV for the potential index u = 1.1. The CPPj,=o.5 
predictions based on V-W with radiative correction is also found to be in good agreement with 
the expected values while in all other choices of i' over estimates the decay width. It indicates 
the importance of the computation from of the decay width at finite range of quark- antiquark 
separation. 

In the case of the leptonic decay width of T(IS') state, most of the models do provide the decay 
widths in close agreement with the expected value either using NRQCD formalism or using V-W 
with radiative corrections. 

Here, again the ERHM prediction for both J/ip and T are found to be very close to the respective 
experimental values with the conventional V-W formula only. It is suggests the adequacy of the 
model parameters that provide the spectroscopy as well as the decay properties. 
To summarize, we find that the spectroscopy of cc system (IS to 3S) studied here are in good 
agreement with the respective experimental values in the potential range of 1.1 < < 1.3. However, 
the spectroscopic predictions with potential index z/ = 1.5 for 66 system are found to be in agreement 
with the respective experimental value. The spectroscopic predictions of the 6c system in the 
potential range 1.1 < z/ < 1.3 are found to be in accordance with other model predictions. 
In the case of the di-gamma decay widths of cc system, better agreement occurs for the potential 
index z/ = 0.7 under the NRQCD and conventional V-W formula with radiative correction. However, 
the NRQCD provide the experimental value of the decay width in the potential index range of 
1.3 < u < 1.5 only. For the 66 system, better consistency in the predictions of both the leptonic 
and di-gamma widths are observed around the potential index 0.7 < z/ < 1.1. 

The present study of the decay rates of quarkonia clearly indicates the relative importance of 
QCD related corrections on the phenomenological potential models. The success of potential mod- 
els in the determination of the 5* and P wave masses and decay rates of cc and 6c and 66 systems 
provide future scopes to study various transition rate and excited states of these mesonic systems. 
With the masses and wave functions of the heavy flavour mesons at hand, it would be rather sim- 
ple to compute various transition rates such as El and Ml in these mesons. Such computations 
largely form the future applications of the present study. The decay rates and branching ratios of 
heavy flavour mesonic bound states are important ingredients in our present understanding of QCD. 

The semileptonic decays offer an extremely favorable testing ground for both perturbative 
QCD, radiative corrections and nonperturbative QCD effects such as decay constants, form factors, 
and the best possible estimations of the CKM matrix elements. With the mass parameters of the 
beauty and the charm quark fixed from the study of its spectra, we have successfully computed the 
semi-leptonic decay width of i?c-meson. 

The partial widths obtained here within the spectator model are compared with those ob- 
tained though the Bethe-Salpeter approach ^24j as well as that from a relativistic quark model [52] 
in Table [TH We obtained a higher branching ratio in the 6-decay channel compared to other ap- 
proaches as seen from in Table [HI We get about 64% as the branching fractions of b-quark decay, 
about 33% as that of c-quark decay and about 3% in the annihilation channel. However, the CKM 
mixing matrix elements Vcb and Vcs used as free parameters in all the three models are different but 
lie within the range given in particle data group |2]. The lifetime of 5+ predicted by the present 
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calculation is found to be in good agreement with the experimental values as well as that by the 
Bethe Salpeter method [23]. The predicted values from relativistic model [52] is found to be far 
from the experimental values as well as other theoretical models. 

Another aspect of the present study is that the decay of QQ system occurs at a finite range of its 
seperation provided by the color compton radius. This enable us to understand at least qualitatively 
the importance of various processes that occur at different radial seperation. 

In conclusion, we have studied the importance of the spectroscopic parameters of different 
potential models in the predictions of the low-lying states of cc, cb and bb systems as well as their 
decay properties in the frame work of NRQCD formalism. 
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